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A hybrid model for chaotic front dynamics: From semiconductors to water tanks 
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We present a general method for studying front propagation in nonlinear systems with a global 
constraint in the language of hybrid tank models. The method is illustrated in the case of semicon- 
ductor superlattices, where the dynamics of the electron accumulation and depletion fronts shows 
complex spatio-temporal patterns, including chaos. We show that this behavior may be elegantly 
explained by a tank model, for which analytical results on the emergence of chaos are available. In 
particular, for the case of three tanks the bifurcation scenario is characterized by a modified version 
of the one-dimensional iterated tent-map. 
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Introduction. Moving fronts are the source of complex 
patterns in very different areas of physics Q] , chemistry 
[S ll| and even biology 0, Q . While the propagation of 
single fronts is well understood for many systems, only 
little is known on how front generation and annihilation 
processes or collisions between fronts influence the possi- 
ble bifurcation scenarios. 

For the sake of concreteness we illustrate our method 
with a semiconductor superlattice system,which is well 
studied experimentally and theoretically 0, IS El EH 
ITU and shows a wide range of complex front patterns (cf. 
Fig-ID including chaos 0,EJE3[- Similar front patterns 
occur in many other systems [T3, e.g., spatially contin- 
uous semiconductor models describing impurity impact 
ionization breakdown |l6l or globally coupled heteroge- 
neous catalytic reactions |l7|. Our aim is to point out a 
striking connection of those front patterns with discrete 
hybrid models of n tanks such as they occur, for instance, 
in industrial production processes [l8| , and derive a sim- 
ple description of the front generation and annihilation 
dynamics in terms of an rt-tank model which allows for 
an analytical treatment of the bifurcation scenarios. In 
particular, an unusual chaotic scenario will be character- 
ized by this generic model. We stress that the presented 
methods do not depend on the pecularities of the super- 
lattice model, and are expected to be equally applicable 
to other front systems with global coupling. 

In order to formulate a generic model for front propa- 
gation, we will first develop a front model based on the 
spatio-temporal dynamics in a semiconductor superlat- 
tice. The dynamical variables are the positions of the 
accumulation and depletion fronts which evolve with a 
velocity determined by the number of fronts present in 
the system. The number of fronts changes according to 
simple rules which mimic generation of fronts at the emit- 
ter and the annihilation of fronts. Such models are quite 
general and similar ideas have been sketched for front 
propagation in bulk GaAs 16]. It will be shown that 
this front model is identical to an n-tank switched ar- 
rival system [Tsl Il9| , if one identifies the spacing Xi (see 
Fig. ^) between the leading depletion and the trailing 
accumulation front with the filling heights of the tanks. 
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FIG. 1: (a) Space-time plot of the electron density evolution 
in a semiconductor superlattice. Electron accumulation and 
depletion fronts are indicated by light and dark shading, re- 
spectively. The full microscopic sequential tunneling model of 
Ref. E| for N =100 quantum wells and contact conductivity 
a — 0.5 Q,~ 1 m~ 1 is used for a bias U = 0.95 V. The inset 
shows a simplified view of the pattern, which is adopted in 
the front and the tank model. The dashed lines denote the 
switching times t m . 



The tanks drain at a constant rate \x, while one of the 
tanks is connected to a server and filled at a rate A = fifi 
(Fig where n is the number of tanks which are not 
empty. In the present model the server will switch to a 
different tank, if this tank is empty and the filling level 
of the currently filled tank is larger than x rn i n . In the 
simplest case n = 3, the tank-model transforms into a 
one-dimensional iterated map which is very similar to 
the well known tent-map. We will show that the front 
dynamics in semiconductor superlattices as simulated in 
1 1 11 ] can largely be explained by such a one-dimensional 
map. Inparticular, atypical border collision bifurcations 
occur [2(| ■ The general type of models we consider here 
is called hybrid since they contain continuous as well as 
discrete state variables. 

Modeling the front dynamics. The superlattice dy- 
namics is described by the electron densities n m (t) in the 
quantum wells, labeled by m = 1, 2 ... TV. The current 
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FIG. 2: Scheme of an n-tank switched arrival system with 
minimal filling height x m in- 
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FIG. 3: Current density vs electric field characteristic at 
the emitter barrier (straight line, Ohmic conductivity a = 
0.5 0,~ 1 m~ 1 ) and between two neutral wells (N-shaped 
curve) for the microscopic superlattice model |14|. J c denotes 
the intersection point of the two characteristics. F l and F h 
are the stable low and high field states, respectively, for the 
current density J c . Inset: Positions where accumulation and 
depletion fronts annihilate vs. voltage. 



density J m -> m +i{n m ,n m+ i, F m ) between adjacent wells 
is evaluated by the sequential tunneling model, where 
F m < is the electric field between these two wells 10]. 
A microscopic calculation of J m ^m+i yields an N-shaped 
current density vs. electric field characteristic as depicted 
in Fig. [3J In the electrically neutral regions of the su- 
perlattice, there are two stable states for fixed current 
density J, one at a low field F l (J) and one at a high field 
F h (J) (Fig|3J. Since in the following J is varied only in a 
small range, we make the approximation that both fields 
do not depend on the current density, and set F l (J) = 0, 
F h (J) = F h > 0. 

The emitter boundary condition (at well m = 1) is 
characterized by its intersection point J c with the charac- 
teristics between neutral wells. For small J c propagating 
inhomogeneous field distributions alternating between F l 
and F h occur, see Fig. JQ. For a transition from a low 
field domain F m — to a high field domain F n = —F h 
(n > m) a negatively charged electron accumulation front 



located at the domain boundary, is required by Poisson's 
equation. Let us denote by et^ the distance of its center 
of charge from the emitter. The index i = 1 . . . N a labels 
different accumulation fronts such that a$ < flj+i, where 
N a is the total number of accumulation fronts in the sys- 
tem. Similarly, a transition from a high field domain 
to a low field domain can be attributed to an electron 
depletion front whose center of charge is located at di 
with i = 1 . . . Na, and Nd the number of depletion fronts. 
Since we only allow for alternating accumulation and de- 
pletion fronts, we have Nd — N a 6 {—1,0,1}. 

The applied voltage U between emitter and collector 
defines the partial length Lh — U/ F h of the superlattice 
which is at high field. This imposes a global constraint 
on the front positions by 



jv,, 



i=i 



mod L. 



(1) 



Here L denotes the total length of the superlattice. The 
expression mod L in means that L has to be added if 
a-N a > dN d such that Lh <E [0, L]. 

The front velocities are given functions of the total cur- 
rent di — v a (J), di = Vd(J), as shown in Refs. |2ll |22|. 
Differentiating eq. QJ yields v d (J)/v a (J) = N a /N d . This 
determines the current J(N a /Nd), which is a monotoni- 
cally increasing function, as v d (J) and v a (J) are increas- 
ing and decreasing, respectively |14| . 



Rescaling the time 



such that 



v d = 2 gives 

2N d 
N a + N d 



2N n 



N a + Nd 



(2) 



Front injection at the emitter is governed by J c . For 
J < J c the region close to the emitter is pinned at a 
low field 01 • Hence an accumulation front is injected, 
if the preceding front is a depletion front. Due to the 
finite width and buildup time of the fronts new fronts 
can not be injected arbitrarily fast. Phenomenologically 
we therefore introduce a distance parameter ph, and as- 
sume that front generation is suppressed while d\ < ph, 
where d\ is the position of the first depletion front. In 
the same fashion, a depletion front is injected if J rises 
above J c and pi < a\, where pi is a parameter describing 
the minimum front distance for depletion front injection. 
As the current J is a monotonic function of N a /N d , the 
condition J ^ J c can be transformed into N a /Nd ^ r c 
with a new critical parameter r c . 

The processes which reduce the number of fronts are 
collisions of two fronts of opposite polarity and the run- 
ning out of fronts at the collector. Both processes affect 
N a /Nd and potentially trigger generation of a new front 
at the emitter. 

Let us summarize the rules for front dynamics: 

Fl The positions of the accumulation fronts for 
i = 1 . . . N a and depletion fronts di for i — 1 . . . Nd 
evolve according to di = v a and di — Vd with the 
velocities until one of the following rules applies. 
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FII If N a /Nd < r c and ph < d± < a\ then increase N a 
by one, re-index a.j — > a^+i for all i and set ai = 
(injection of accumulation front). 

Fill If N a /Nd > r c and pi < a\ < d\ then increase N d 
by one, re-index di — > <fj+i for all i and set d\ = 
(injection of depletion front). 

FIV If Oj' = dj' for any i',f then decrease iV a and N d 
by one, re-index a,+i — > ai for i > i' and — > dj 
for j > j' (annihilation of front pair). 

FV If a,N a > L decrease N a by one (accumulation front 
hits collector). 

FVI If dN d > L decrease Nd by one (depletion front hits 
collector). 

Note that the voltage parameter Lh only enters in the 
initial condition for the front positions (see Eq. QJ). 

In the following we will restrict ourselves to the case 
r c < 1 and pi = 0. Rule FII does not apply as long as 
N a > r c /(l — r c ) and N a can then only decrease. Conse- 
quently r c imposes the following limits on the number of 
fronts: 

N d <n; N a < n - 1; (3) 

where n is the largest integer less than 1/(1 — r c ) + 1. 
The injection of an accumulation front is immediately 
followed by the injection of a depletion front (Fill). This 
detaches a high field domain from the emitter and opens 
a new one. Furthermore we restrict ourselves to the case 
where no fronts reach the collector (i.e. L > nLh). This 
corresponds to the situation in Fig. JTJ and we have iV a = 
N d - 1. 

The n-tank model. Under the restrictions given above 
we may choose the lengths of the high field domains (cf. 
inset of Fig.l) x\ = d\ and Xi = di — a,_x for i — 2 . . . Nd 
as new dynamical variables governed by the rules: 

TI The lengths evolve according to ii = —p, + A<5,i 
with fi = 2/(2N d - 1) and A = N d fi. 

Til If Nd < n and x\ > ph then increase N d by one, 
re-index Xi — > Xi+\ for all i and set x± = 0. 

Till If Xi> = then decrease N d by one, re-index a;,_|_i — » 
Xi for all i > i' . 

The voltage constraint is expressed in the new vari- 
ables by J2k x k = Lh- 

It is very illustrative to interpret this model in the 
language of an n-tank switched arrival system (cf. FigJ2J. 
A domain may be viewed as a tank with fluid content Xi. 
One tank is filled at rate (Nd — while N d — 1 tanks 
are drained at rate [A, and n — Nd tanks are empty and 
waiting for filling, because the currently filled tank has 
not yet reached the minimum filling level x m i n = ph- 

Let us first discuss the limiting case pi — Ph = 0. 
If a tank is empty at time t m (i.e. Xk(t m ) — 0), the 
rules Til and Till apply instantaneously and the system 




FIG. 4: Bifurcation diagram of xi vs. Lh in units of ph- 
Inset: Graph of the return map for the n — 3 tank model Eq. 
©. In the shaded region the map is not defined. 

starts to fill this tank (aq(i+) = 0). In this case no 
tank has to wait for filling. Apart from the re-indexing 
procedures we obtain a switched arrival system |l£j |. As 
shown in [2j| the system is chaotic for all n > 2 and has 
a constant invariant probability measure. By sampling 
the dynamics at the switching times t m we obtain the 
Poincare map P : x(t m ) i— ► x(t m +i): 

x{t m+ x) = x(t+ n ) + SAt m (4) 

where x = (A — ji, —ji, . . . , — u) is a constant vector and 
At m = t m+ i — tm = mizxj^\(xj(t m )/ n). Generally this is 
an n dimensional mapping which acts due to the normal- 
ization Xi = Lh on the n — 2 dimensional boundary of 
a regular n-simplex embedded in R™. For n = 3 we can 
reduce this map to a one-dimensional map of the interval 
[0, Lh] onto itself: From the rules we obtain £i(f+) = 0, 
x 2{t m ) = %i(t m ), and ar 3 (t+) = L h -Xi(t m ). If Xi(t m ) < 
Lh/2, the tank 2 empties after At m = Xx(t m )/pb and we 
find xi(t m+ i) = 2xi(t m ). Otherwise, if xi(t m ) > Lh/2, 
the tank 3 empties after At m — (Lh — xi(t m )) / /j, and 
we find xi(t m +i) — 2(Lh — x\(t m )). Thus the dynam- 
ics is that of the well known tent map: x\(t m+ i) = 
L h -2\x 1 (t m )~L h /2\. 

Now we consider the case n — 3 with a finite minimum 
filling level ph > 0. Since by Til x\ must be larger than 
Ph for switching, the Poincare map is now only defined 
on the interval [ph,Lh]- For X\ £ [ph, Lh —ph/2] the map 
remains unchanged from the tent map discussed before. 
If the switch to an empty tank occurs at x\(t m ) > Lh — 
Ph/2, the third tank will be empty before the newly filled 
tank has reached the minimum filling level ph- Switching 
is therefore suppressed until xi(t m+ i) = ph- Thus for 
n = 3 we obtain (cf. inset of FigQJh 

a?i(Wi) = Max{(i h - 2\ Xl (t m ) - L h /2\), Ph } (5) 

The bifurcation structure of this map with respect to Lh 
is unusual (Figffl but has similarities to the flat-topped 
tent map |24i l2q . Due to the horizontal segments the 
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map cannot show chaotic behaviour although it under- 
goes several bifurcation cascades. A trajectory on the 
tent map segments would explore the attractor of the 
tent map and will therefore either fall on one of the flat 
segments or stay at an unstable periodic orbit contained 
in the remaining part of the tent map. Once arrived at 
one of the flat segments the trajectory continues on a 
stable periodic orbit. By analyzing the fc-th return map 
of the modified tent map we find that windows of period 
k are found inside the intervals L h /p h £ { 2 A jZ\ > §7^5 ) f° r 
j e h C {j E N\j < 2 k ~ 2 - 1}. Thus higher and higher- 
stable periodic orbits are arising in smaller and smaller 
intervals. 

Discussion. In order to compare the bifurcation dia- 
gram of the modified tent map (FigQJ with the bifurca- 
tion diagram of the superlattice model (inset of Fig. |2J , 
we note that the voltage U and the position of front an- 
nihilation in the superlattice model correspond to the 
parameters Lh and xi(t m ) in the map, respectively. We 
can therefore identify the periodic windows found in the 
Poincare map with the periodic pattern observed in the 
front model. We note that the bifurcation scenarios 
agree very well up to voltages of about 1.1V, where a 
period-3 window occurs (inset of Fig. |3J). In particular, 
the remarkable spider-like structure at 0.9V is found for 
Lh/Ph = 3 in the bifurcation diagram of the map. The 
reason for this spider is that the contracting flat region is 
mapped onto the unstable fixed point of the map. Fur- 
ther prominent features in both bifurcation diagrams are 
the horizontal lines given by the first, and higher iterated 
images of the minimal filling height. Thus the character 



of the bifurcation scenario of complex front dynamics is 
explained by considering a one-dimensional iterated map 
with only one free parameter Lh/ph- 

The minimum filling height corresponds to the mini- 
mum domain width. For a real superlattice system this 
quantity will slightly depend on the history of the dynam- 
ical evolution. Therefore the slope of the flat segments is 
not exactly zero but attains a small finite value a. In this 
case a period-fc orbit is unstable if \2 k a\ > 1. Close to 
the spider-like points where periodic orbits with large pe- 
riod k occur, we can then find a finite chaotic regime, as 
observed in our simulations with the microscopic model. 

Conclusions. We have introduced a hybrid tank 
model as a simplified model for front dynamics in non- 
linear systems with a global constraint. This simplified 
model provides an explanation for the complex accumu- 
lation and depletion front dynamics in semiconductor su- 
perlattices. In this unexpected analogy, the tank filling 
levels correspond to the widths of the high field domains 
in the superlattice, while the position of the server deter- 
mines which high field domain is connected to the emit- 
ter. Although this analysis has started from a superlat- 
tice model, the presented methods are very general and 
are expected to work successfully for other systems which 
provide a global coupling similar to Eq. 
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